RANDOM REPEATED QUANTUM INTERACTIONS 
AND RANDOM INVARIANT STATES 



ION NECHITA AND CLEMENT PELLEGRINI 

Abstract. We consider a generalized model of repeated quantum interactions, where a sys- 
tem Ti is interacting in a random way with a sequence of independent quantum systems 
IC„,n ^ 1. Two types of randomness are studied in detail. One is provided by considering 
Haar-distributed unitaries to describe each interaction between Ti. and IC„. The other in- 
volves random quantum states describing each copy /C„. In the limit of a large number of 
interactions, we present convergence results for the asymptotic state of 7i. This is achieved by 
studying spectral properties of (random) quantum channels which guarantee the existence of 
unique invariant states. Finally this allows to introduce a new physically motivated ensemble 
of random density matrices called the asymptotic induced ensemble. 



1. Introduction 

Initially introduced in [2] as a discrete approximation of Langevin dynamics, the model 
of repeated quantum interactions has found since many applications (quantum trajectories, 
stochastic control, etc.). In this work we generalize this model by allowing random interactions 
at each time step. Our main focus is the long-time behavior of the reduced dynamics. 

Our viewpoint is that of Quantum Open Systems, where a "small" system is in interaction 
with an inaccessible environment (or an auxiliary system). We are interested in the reduced 
dynamics of the small system, which is described by the action of quantum channels. When 
repeating such interactions, under some mild conditions on the spectrum of the quantum 
channel, we show that the successive states of the small system converge to the invariant 
density matrix of the channel. 

These considerations motivated us to consider random invariant states, and we introduce 
a new probability measure on the set of density matrices. There exists extensive literature 
[HI dSl [3] on what is a "typical" density matrix. There are two general categories of 
such probability measures on Ai^'~^{C): measures that come from metrics with statistical 
significance and the so-called "induced measures", where density matrices are obtained as 
partial traces of larger, random pure states. Our construction from Section [J] falls into the 
second category, since our model involves an open system in interaction with a chain of 
"auxiliary" systems. 

Next, we introduce two models of random quantum channels. In the first model, we allow 
for the states of the auxiliary system to be random. In the second one, the unitary matrices 
acting on the coupled system are assumed random, distributed along the Haar invariant 
probability on the unitary group, and independent between different interactions. Since the 
(random) state of the system fluctuates, almost sure convergence does not hold, and we state 
results in the ergodic sense. 
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The article is structured as follows. The Section [2] is devoted to presenting the model of 
quantum repeated interactions and its description via quantum channels. Section [3] contains 
some general facts about the spectra of completely positive maps, as well as some related 
tools from matrix analysis. Next, in Section 5] we study our first model, where the interaction 
unitary is a fixed, deterministic matrix. We prove that, under some assumptions on the 
spectrum of the quantum channel, the state of the system converges to the invariant state 
of the channel. It is at this time that we introduce the new ensemble of random density 
matrices, by transporting the unitary Haar measure via the application which maps a channel 
to its invariant state. The final two sections are devoted to introducing two models of random 
quantum channels, one where the interaction unitary is constant and the auxiliary states are 
i.i.d. density matrices (Sec. [5|) and another where the interaction unitaries are independent 
and Haar distributed (Sec. [6]). 

We introduce now some notation and recall some basic facts and terminology from quantum 
information theory. We write A^^'^(C) for the set of self-adjoint d x d complex matrices and 
MI'^{C) for the set of density matrices (or states), M^'^iC) = {p e Mf{C) \ p ^ 0,Tr[p] = 
1}. Since our main focus is quantum information, all Hilbert spaces in this article are complex 
and finite dimensional. Scalar products are assumed linear in the second coordinate and, for 
two vectors x & Tl,y G }C we denote by \x){y\ € B{}C,H) the map 



An unit vector x G W ~ is called a pure state and it is assimilated often with the orthogonal 
projection on Cx, \x){x\. Finally, for a matrix A G 0(7i (g) /C) ~ 13{TC) B{1C), we define its 
partial trace with respect to K, as the unique element B = Tt)q[A\ G B{T-L) which verifies 



We shall also extensively use the Haar (or uniform) measure 1)^ on the unitary group lA{d); 
it is the unique probability measure which is invariant by left and right multiplication by 
unitary elements: 



In this introductory section we give a description of the physical model we shall use in the 
rest of the paper: repeated quantum interactions. The setting, a system interacting repeatedly 
with "independent" copies of an environment, was introduced by S. Attal and Y. Pautrat in 
[2] where it was shown that in the continuous limit (when the time between interactions 
approaches zero), the dynamics is governed by a quantum stochastic differential equation. A 
different model, where after each interaction an indirect quantum measurement of the system 
is performed, was considered by the second named author in [211 [22] and shown to converge 
in the limit to the so-called stochastic Schrodinger equations. Here, we are concerned only 
with the discrete setting and with the limit of a large number of interactions. The study of 
random quantum trajectories is postponed to a later paper. 

Consider a quantum system <S described by a complex Hilbert space state 7i. In realistic 
physical models, S is usually a quantum system with relatively few degrees of freedom and it 
represents the object of interest of our study; we shall refer to it as the small system. Consider 
also another quantum system £ which interacts with the initial small system iS. We shall call 
£ the environment and we denote by /C its Hilbert state space. In this work we consider finite 
dimensional spaces 7^ ~ and /C ~ . 



x){y\{z) = {y,z) ■ X, yzeJC. 



Tv[BX] = Tv[A{X (g> lie)], VX G B{n). 



VF, W G U{d), V/ : U{d) C Borel, 




2. The repeated quantum interaction model 
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We shall eventually be interested in repeated interactions between S and independent copies 
of £, but let us start with the easier task of describing a single interaction between the "small" 
system and the environment. Assume that the initial state of the system is a product state 
a = p CS) P, where p and /3 are the respective states of the small system and the environment. 
The coupled system undergoes an unitary evolution U and U{p (3)U* is the global state 
after the interaction. The unitary operator U comes from a Hamiltonian 

Htot = Hs®l + l ®H£ + Hint, 

where the operators Hg and are the free Hamiltonians of the systems S and £ respectively 
and Hint represents the interaction Hamiltonian. We shall be interested in the situation where 
Hint 7^ 0, otherwise there is no coupling and the system and the environment undergo separate 
dynamics. In this general case, the evolution unitary operator U is given by 

U = e-"^'"', 

where r > is the interaction time. Hence, the state of the coupled system S + £ after one 
interaction is given by 

a' = U{p®())U*. 

Since one is interested only in the dynamics of the "small" system 5, after taking the partial 
trace we obtain the final state of 5, 

(1) p' = T^K.[u{p®m*]- 

We now move on to describe successive interactions between S and a chain of independent 
copies of <f^. In order to do this, consider the countable tensor product 

oo 

^tot = (^^n, 
n=l 

where /C„ is the n-th copy of the environment (/C„ ~ /C ~ C^' ) . This setting can be interpreted 
in two different ways: globally, as an evolution on infinite dimensional countable tensor prod- 
uct Ti ® fCtot, or by discarding the environment, as a discrete evolution on B{Ti.) = A^d(C). 
Since we are interested only in the evolution of the "small" system, the latter approach is the 
better choice. Prom Eq. ([T]), we obtain the recurrence relation 

(2) Pn = Tr,c[Un{Pn-10Pn)K], 

where pn-i, Pn G -^d'^lC) are the successive states of the system S at times n — 1 and n, and 
Un and Pn are the interaction unitary and respectively the state of the auxiliary system £ for 
the n-th interaction. Note that at this stage we work in a general setting, without making 
any assumptions on the sequences {Un)n and {Pn)n- 

We introduce now a more parsimonious description of repeated quantum interactions, via 
quantum channels. Recall that a linear map ^ : A4d{C) — > A^rf(C) is called /c-positive if 
the extended map $ (g) Ifc : Md{C) (g) MkiC) MdiC) Mk{C) is positive. $ is called 
completely positive if it is /c-positive for all A; ^ 1 (in fact k = d suffices) and trace preserving 
if Tr[<I>(X)] = Tr[X] for all X G A4rf(C). By definition, a quantum channel is a trace- 
preserving, completely positive linear map. The next proposition gives two very important 
characterizations of quantum channels. 

Proposition 2.1 (Stinespring-Kraus). A linear map <I> : A4^(C) A4rf(C) is a quantum 
channel if and only if one of the following two equivalent conditions holds. 
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(a) (Stinespring dilation) There exists a finite dimensional Hilbert space IC = , a 
density matrix [5 G AA^^iC) and an unitary operation U G lA{dd') such that 

^X) = Trc [U{X ® 0)U*] , G Md{C). 

(b) (Kraus decomposition) There exists an integer k and matrices Li, . . . , G MdiC) 
such that 

k 

(3) ^x) = Y,LiXL*, yxeMdiC) 

and 

k 

= Id . 

i=l 

Remark 2.2. It can be shown that the dimension of the ancilla space in the Stinespring dilation 
theorem can be chosen = d'^ and /? can be chosen to be a rank one projector. A similar 
result holds for the number of Kraus operators: one can always find a decomposition with 
k = operators. The Choi rank of a quantum channel <I> is the least positive integer k such 
that ^ admits a Kraus decomposition ([3]) with k operators Lj. 

We see now that Eq. ([2|) can be re-written as 

Pn = $^"'''"(pn-l), 

where is the quantum channel 

/J ^Tr/c[f/(p® /?)[/*]. 
After n such interactions, the state of the system becomes 

(4) Pn = O 0---0 ^'^I'^l/J- 

Let us now consider a fixed channel <I> = and show that the Stinesping and Kraus 
form of ^ are connected in a simple fashion. To this end, start with the Stinespring form of 
$ and pick some orthonormal bases {ei}f^i and {fj}j=i of respectively TL = and K. = 
such that the state of the environment /3 diagonalizes: 

d' 

Next, endow the product space n®IC = C^' with the basis 

(5) {ei (g) /i, 62 /i, . . . , e„ (g) /i, ei (g) /2, . . . , Cn O /2, • • • , e„ fk}- 

This particular ordering of the tensor product basis was preferred in order to have a simple 
expression for the partial trace operation with respect to the environment /C. Indeed, if 
a matrix A G Mdd'i'C) is written in this basis and viewed as a d' x d' matrix of blocks 
Aij G Md{C): 

/All Ai2 ••• Aid'\ 

A2I A22 ■ ■ ■ A2d' 



A 



X^d'i Ad'2 ■ ■ ■ Ad'd' J 
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then the computation of the partial trace with respect to /C = reads 
Tt,c[A]=Ttk 



«J=1 



d' 

A,j ■ {fj, /,) = All + A22 + • ■ ■ + Ad'd' 



In other words, the partial trace of A over the environment /C is simply the trace of the 
block-matrix, that is the sum of the diagonal blocks of A. We apply now these ideas to the 
Stinespring form of a quantum channel, ^{X) = Tr k^[U {X P)U*]. Written as a block matrix 
in the basis defined in Eq. ([5]), the matrix X (8) /3 is diagonal, with diagonal blocks given by 
bjX € A4(i{C). Writing U £ U{dd') in the same fashion and taking the partial trace, we 
obtain 

d> d' 

(6) $(X) = Tt^[U{X ® /?)[/*] = Y bjU.jXU:^ = Y {v^U,,)X{y^U,,r, 

i,j=l i,j=l 

where Uij G A^rf(C) are the blocks of the interaction unitary U. One recognizes a Kraus 
decomposition for where the Kraus elements are rescaled versions of the blocks of the 
Stinespring matrix U. Moreover, if /? is a rank one projector then all the 6j's are zero except 
one, hence the Kraus decomposition we obtained has d' elements. 

3. Spectral properties of quantum channels 

Since we shall be interested in repeated applications of quantum channels, it is natural that 
spectral properties of these maps should play an important role in what follows. One should 
note that most results of this section can be generalized to infinite dimensional Hilbert spaces. 

The next lemma gathers some basic facts about quantum channels. Since quantum channels 
preserve the compact convex set of density matrices AI^'^(C), the first affirmation follows from 
the fixed point theorem of Markov-Kakutani p!0]- The second and the third assertions are 
trivial (see [23] for further results on norms of quantum channels), and the last one is a 
consequence of 2-positivity. 

Lemma 3.1. Let $ : A4d{C) 7V4rf(C) a quantum channel. Then 

(a) $ has at least one invariant element, which is a density matrix; 

(b) ^ has trace operator norm of 1; 

(c) $ has spectral radius of 1; 

(d) $ satisfies the Schwarz inequality 

VX G Md{C), $(X)*$(X) ^ ||«>(1)|| <^{X*X). 

If one looks at a channel $ as an operator in the Hilbert space 7Wrf(C) endowed with the 
Hilbert-Schmidt scalar product, then one can introduce ^, the dual map of It is defined 
by the relation 

Tr[X$(y)] = Tt[^{X)Y], VX, Y G Md{C). 

From Kraus decomposition ^{X) = ^ LiXL*, one can obtain a Kraus decomposition for the 
dual channel, ^{X) = ^ L*XLi. Note that the trace preserving condition for ^ L*Li = I 
reads now ^'(I) = I. Hence, the dual of a quantum channel is a unital (not necessarily trace- 
preserving) completely positive linear map. Using this idea, one can see that the partial trace 
operation Tr^; : Mdd'i'C) 7Wd(C) is the dual of the tensoring operation Sjc ■ M.d{C) 
Mdd'iC), S{X) = X0ld'. 

We now introduce some particular classes of positive maps which are known to have inter- 
esting spectral properties. 
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Definition 3.2. Let $ : A4d{C) A^d(C) be a positive linear map. $ is called strictly 
positive (or positivity improving) if ^{X) > for all X ^ 0. $ is called irreducible if there is 
no projector P such that $(P) ^ XP for some A > 0. 

Example 3.3. Let U G Z//((i) be a fixed unitary and consider the channel $ : A4d{C) A^rf(C), 
$(X) = UXU*. It is easy to check that the spectrum of $ is the set 

{A1A2 I Ai, A2 eigenvalues of U}. 

Since $ maps pure states (i.e. rank-one projectors) to pure states, it neither irreducible, nor 
strictly positive. 

Obviously, a strictly positive map is irreducible. In fact, the following characterization of 
irreducibility is known 

Proposition 3.4. A positive linear map $ : Mdi^) A4d{C) is irreducible if and only if 
the map (1 + ^)'^^^ is strictly positive. 

Irreducible unital maps which satisfy the Schwarz inequality have very nice peripheral 
spectra. The proof of the following important result can be found in one of [111 [T3l I15j . in 
more general settings. 

Theorem 3.5. If^ is a unital, irreducible map on A4d{C) which satisfies the Schwarz inequal- 
ity, then the set of peripheral (i.e. modulus one) eigenvalues is a (possibly trivial) subgroup 
of the unit circle T. Moreover, every peripheral eigenvalue is simple and the corresponding 
eigenspaces are spanned by unitary elements of Aid{C). 

Irreducible (and, in particular, strictly positive) quantum channels have desirable spectral 
properties, hence the interest one has for these classes of maps. As we shall see in Section 
m irreducible maps are in certain sense generic. On the other hand, the strict positivity 
condition is rather restrictive and not suitable for the considerations on this work. Next, we 
develop these ideas, giving criteria for irreducibility and for strict positivity. 

Let us start by analyzing strict positivity. Subspaces of product spaces with 
high entanglement have received recently great attention. In this direction, applications to 
the additivity conjecture \16\ [T7] are the most notable ones. The results in these papers, 
which rely on probability theory techniques deal with von Neumann entropy. When one looks 
at the rank, projective algebraic geometry comes into play. Indeed, possible states of the 
coupled system are modeled by the projective space space contains the product 

states, P'^"! P'^ ~^ as a subset called the Segre variety. The following lemma, a textbook 
result in algebraic geometry, is obtained by computing the dimension of the Segre variety (see 
13 [201 [26]). 

Lemma 3.6. The maximum dimension of a subspace S C <^ C^' which does not contain 
any non-zero product elements x ® y is {d — l)(d' — 1). 

As a rather simple consequence of this lemma, we obtain a necessary condition for strict 
positivity. 

Proposition 3.7. Let <I> : A^rf(C) M.di'^) strictly positive quantum map. Then the 
Choi rank of <I> is at least 2d — 1. 

Proof. Let ^{X) = LiXL* be a minimal Kraus decomposition of a strictly positive 

channel <I>. For all a: 7^ 0, $([x)(x|) has full rank, and thus, for all non-zero y £ C^, 

k 

TTmx){x\)\y){y\]=Y,\{y,Lix)\^>^. 

i=l 
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Hence, for all non-zero x,y £ C^, there exist an i such that {y,Lix) = Tr[Lj[x)(y|] ^ 0, 
or, in other words, L* is not orthogonal to with respect to the Hilbert-Schmidt scalar 

product. Consider now the space S = C]i=im)'^ C MdiC). Obviously, S does not contain 
any rank one matrices \x){y\. Under the usual isomorphism (C^)* ~ A4d{C), product 
vectors x®y are identified with rank one matrices \x){y\, so, by the Lemma 13.61 we get 
dimS ^ {d — 1)^ = — {2d — 1). Since S is the intersection of k subspaces of dimension 
- 1, we get - /c ^ dim S d'^ - {2d - 1) which imphes k^2d-l. □ 

We now turn to irreducible quantum maps and state some results which will be useful later, 
when showing that irreducibility is generic for a specific model of random quantum channels. 

The following result of [13] gives necessary and sufficient conditions for a map written in 
the Kraus form to be irreducible. We denote by Lat(T) the lattice of invariant subspaces of 
an operator T G Md{C). 

Proposition 3.8. Consider the map : Md{C) Md{C) defined by ^{x) = Yli=i LiXL*, 
with Li € M.d{'C-), i = 1, . . . , k. Then ^ is irreducible if and only if C\i=i Lat(Li) is trivial. 

Of course, quantum channels of Choi rank one (i.e. unitary conjugations, see also Example 
13. 3p . ^{X) = LXL*, with L*L = I cannot be irreducible, since they leave invariant eigen- 
projectors of L. When looking at channels with Choi rank at least two, an useful criterion 
for deciding whether 0^=1 Lat(Lj) is trivial or not is given by the following two results. The 
first proposition gives necessary and sufficient conditions for two matrices A and B to share a 
common eigenvector, and the second one generalizes this idea to arbitrary common subspaces. 

Proposition 3.9 (The Shemesh criterion, |24j). Two matrices A, B £ ^Ad{C) have a common 
eigenvector if and only if 

d-l 

f] keT[A\B^]^{0}, 

or, equivalently, iff 

d-l 

det ^ • S-''] = 0. 

In order to move on from common eigenvectors to common invariant subspaces, we consider 
antisymmetric tensor powers (or wedge powers) of matrices (see [4], Ch. I). Given A € Md{C) 
and and an integer 1 ^ fc ^ n, the A;-th wedge power of A, denoted by A'^^ , is defined as 
the restriction of A®^ to the antisymmetric tensor product (C"')^'^. More precisely, A^^ is a 
n X n matrix, where n = . Its matrix elements are indexed by couples {a, (3) of strictly 
increasing sequences of size k from {1, . . . , d}: 

(a'^'') =detyl[a|/3], 

where is the submatrix of A with rows indexed by a and columns indexed by [3. The 

next result of |14j is an easy consequence of the fact that if Ai,...,Afc are eigenvalues of 
A with linear independent vectors wi, . . . then A1A2 • • • is an eigenvalue of A^^ with 
corresponding eigenvector vi /\ ■ ■ ■ /\Vk- 

Proposition 3.10 (Generalized Shemesh criterion, [13]). Let A,B £ Md{'C-) be two complex 
matrices. If A and B have a common invariant subspace of dimension k (for 1 ^ k ^ d — 1), 
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then their k-th wedge powers have a common eigenvector, and hence (we put n 

n-l 

n ker[{A^')\iB^y]^{0}, 

or, equivalently. 



it)) 



n-l 

det [{^''''y, (B'^'^y] 



0. 



Remark 3.11. The preceding conditions turn out to be sufficient under more stringent as- 
sumptions on the matrices A and B (see [14J for further details). 

The main point of the two preceding results is that there exists an universal polynomial 
P £ . . . ,^4^2] with the property that whenever two matrices A = {aij) and B = {h^i) 

have a non-trivial common invariant subspace, P(Reajj,Imaij,Re6fc;,Im5fci) = 0. This fact 
(together with Proposition 13. 8p will be useful later in this work, when we shall show that a 
generic class of quantum maps are irreducible. 

4. Non-random repeated interactions and a new model of random density 

matrices 

In this section we consider repeated interactions with a fixed unitary matrix U (Vn, C/„ = U) 
and fixed state of the environment (3 {Vn, /3„ = /?). By the results of the previous section, the 
recurrence relation which governs the discrete, deterministic dynamics is 

Pn+l = ^(/9„) = Tlic [U{Pn ® fi)U*] ■ 

Iterating this formula, one obtains the state of the system after n interactions: 

Pn = ^"(po), 

where pQ was the initial state of the system. There is one obvious situation in which the 
asymptotic properties of the sequence {pn)n can be established. Indeed, from Lemma [3.H one 
knows that all quantum channels have eigenvalue 1 and that all other eigenvalues have module 
less than 1. Let C be the set of all quantum channels that have 1 as a simple eigenvalue and 
all other eigenvalues are contained in the open unit disc. Since 1 is a simple eigenvalue, $ has 
an unique fixed point which is (by Lemma IXTI) a density matrix p^ G A^^'^(C). Using the 
Jordan form of <I>, one can show the following result ([25]). 

Proposition 4.1. Let ^ £ C be a fixed quantum channel. Then, for all pQ G A^^'^(C), 

lim ^>"(po) = Poo, 

n^oo 

where poo is the unique invariant state of^. 

The importance of the peripheral spectrum of a quantum channel is illustrated in the 
following example. 



Example 4.2. Consider the following channel ^ : A^2(C) 

where the Pauli matrices are given by 
(To = I 



A^2(C) 



"1 


0' 




"0 


1" 




"0 -i 




"1 " 





1 


, en = 


1 





, 0-2 = 


i 


, 0-3 = 


-1 
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Direct computation shows that <I>(I) = I, <^(o"2) = — (T2 and <&(o"i) = ^{(13) = 0. Hence, the 
peripheral spectrum of $ has 2 simple eigenvalues, 1 and —1. However, for pQ = 1/2(1 +£12) G 
mI'^{C), one has 

^"(po) = ^(i+(-ir^2), 

which does not converge in the limit n ^ 00. Hence, the simplicity of the eigenvalue 1 does not 
suffice to have convergence to the invariant state. Note also that the channel <I> is irreducible, 
since ai and do not have any common non-trivial invariant subspaces. 

Let us now show that the class C of quantum channels which have 1 as an unique peripheral 
eigenvalue is generic in a certain sense. To this end, we shall introduce a model of random 
quantum channel, based on the Stinespring decomposition. To start, fix the dimension of 
the environment d' and a state f3 € Ai^,^{C). Next, consider an unitary random matrix U 
distributed along the (uniform) Haar measure i)dd' on lA{dd'). To the state /? and the evolution 
operator [/, we associate the quantum channel ^u,p- In this way, we define a model of random 
quantum channels by considering the image measure of the Haar distribution \]dd' on the set of 
quantum channels. In the recent preprint [2^, the authors study a similar model of random 
quantum channels, focusing on the spectral properties of the random matrix defining the 
channel. 

More precisely, we claim that if the state of the environment (3 is fixed and the interaction 
unitary U € lA{dd') is chosen randomly with the uniform Haar distribution \]dd' ■, then, with 
probability one, the channel admits 1 as the unique eigenvalue on the unit circle. Here we 
need another fact from algebraic geometry, summarized in the following lemma (for a similar 
result, one should have a look at Proposition 2.6 of [1]). 

Lemma 4.3. Given a polynomial P G M[Xi, . . . , ^2^2], the set 

Z = {U = (uij) G U{d) I P{Reuij,lmuij) = 0} 
is either equal to the whole set lA{d) or it has Haar measure 0. 

Proof. We start by noticing that the real algebraic set U{d) is irreducible. This follows from the 
connectedness oiU{d) (in the usual topology) and from the fact that irreducible components of 
a linear algebraic group are disjoint ([18j, 7.3). The set Z is the intersection of the irreducible 
variety U{d) with the variety V of zeros of the polynomial P. liU{d) C V, then Z = U{d); 
otherwise, the dimension of Z is strictly smaller than d^, the real dimension ofl4{d). Since 
the Haar measure is just the integration of an invariant differential form, it has a density in 
local coordinates (p3], Ch. 5) and hence ijdiZ) = in this case. □ 

Theorem 4.4. Let P be a fixed density matrix of size d' . If U is a random unitary matrix 
distributed along the Haar invariant probability \)dd' on U{dd'), then (^^'1^ G C almost surely. 

Proof. The proof goes in two steps. First, we show that ^'^'^ is almost surely irreducible and 
then we conclude by a simple probabilistic argument. 

Let us start by applying Lemma 14.31 to show that a random quantum channel is almost 
surely irreducible. To this end, using Eq. ([6]), we obtain a set of Kraus operators for <I)^'^ 
which are sub-matrices of [/ G U{dd'). Consider two such Kraus operators A,B A4d{C) 
(choose j such that bj 7^ and take A = Uij, B = U2j). Using Proposition 13.81 to show 
irreducibility it suffices to see that A and B do not have a non-trivial common invariant 
subspace. Let 1 ^ k ^ d — 1 he the dimension of a potentially invariant common subspace 
of A and B. By the criterion in Proposition 13. lOl there exists a polynomial Pk in the entries 
of A and B (and thus in the entries of U) such that if Pk{U) is non-zero, then A and B do 
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not share a A;-dimensional invariant space. Note that Pk can not be identicahy zero: for two 
small enough matrices A, B without common invariant subspaces, one can build a unitary 
matrix U such that A = Uij, B = U2j- By the Lemma 14.31 f)(i(i/-almost all unitary matrices 
U give Kraus operators A and B that do not have any fc-dimensional invariant subspaces 
in common. Since the intersection of finitely many full measure sets has still measure one, 
almost all quantum channels are irreducible. 

Consider now a random channel which we can assume irreducible. Since the peripheral 
spectrum of an irreducible channel is a multiplicative subgroup of the unit circle T, it suffices 
to show that for all element A of the finite set G T|31 ^ n ^ d? s.t. = 1}\{1}, with Haar 
probability one, A is not an eigenvalue of <I)^'^. We use the same trick as earlier. Consider 
such a complex number A and introduce the polynomial Q\{U) = det[<I>^'^ — AI(£;^/')2], where 
^>'^'^ is seen as a matrix ^'^'^ G 7W(ci(^/)2(C). Since A 7^ 1 and the identity channel has 
only unit eigenvalues, Qx{U) cannot be identically zero, and the conclusion follows. □ 

Remark 4.5. The main difficulty in the proof of the preceding result comes from the fact that 
the matrices A and B are "correlated": two blocks of an unitary matrix must satisfy norm 
and (maybe) orthogonality relations. Hence the need to use sophisticated geometric algebra 
techniques. Proving that two independent random Gaussian (or unitary) matrices do not 
share non-trivial invariant subspaces is much simpler and does not require the use of such 
techniques. 

We now move on and apply the previous results to constructing a new family of probability 
distributions on the set of density matrices. The main idea is to assign, whenever possible, 
to a random unitary U € lA{dd') its unique invariant density matrix p^o. In this way, the 
Haar measure \]dd' on the unitary group lA{dd') is transported to the set of density matrices 
mY{C). 

Let us now make this construction more precise. The new family of probability measures 
shall be indexed by an integer d' ^ 1 (the dimension of the auxiliary system) and by a non- 
increasing probability vector b = {bi, . . . ,bd') € : 61 ^ 62 ^ ... ^ bd' ^ and Y2i = 1 
(these are the eigenvalues of the state of the auxiliary system). For such a couple (d',6) 
consider a density matrix P € A4^;^(C) with eigenvalue vector b (the eigenvectors of /3 do 
not matter, see Lemma 14. 6p . As it follows from Proposition 14.41 for almost all unitaries 
U G lA{dd'), the channel ^'^'^ satisfies the hypotheses of Proposition 14.11 Hence, for almost 
all U and for all density matrices po G Ad^~^{C), Ymin^oo{^^'^)^ Po = Poo, where poo is the 
unique invariant state of <I)^'^. We have defined almost everywhere an application 

U{dd')^MY{C) 

U ^ Poo- 

We denote by vi, the image measure of the Haar probability \]dd' on U{dd') by the previous 
application (notice that we dropped the integer parameter d', since this is the dimension of 
the vector b). We call Vb the asymptotic induced measure on the set of density matrices. 

We now motivate the term "asymptotic induced" in the previous definition by showing 
how the measures Vb relate to the induced random density matrices considered in [28^ I19j . 
Let us recall here how these measures are constructed and how one can sample from this 
distribution. The physical motivation behind the induced measures comes from the following 
setup. Assume that a system S is coupled to an environment £ and that the whole is in a 
pure state ip ^ 7i ® KL. If one has no a priori knowledge about the state V'j then it is natural 
to assume that ip is a. random uniform element on the unit sphere of the product space 7i®lC. 
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The distribution of the partial trace over the environment 

is called the induced measure and it is denoted by ^.d' (the parameter d' = dim/C is the 
dimension of the environment). We refer the interested reader to [19j for more information 
on these measures. Since the distribution of an uniform norm-one vector ip is the equal to 
the distribution of UtpQ^ where tpQ is any fixed norm-one vector and C/ is a Haar unitary, /i^/ 
is also the distribution of the matrix 

pi = TiK.[U\i,o)mU*]. 

If one chooses "00 = ei (8) /i, where ei and /i are the first vectors of the canonical basis of 
and respectively , then 

with pq = |ei)(ei| and /3o = j/i)(/i|- Hence, the induced measure p.^' is the distribution of 
the result of one application of a random channel on the constant matrix pq: pi ^ p^' ■ 

On the other hand, after a large number of identical interactions, one gets 



lim 



+ 00 



(Po). 



In this work we have shown that with fj^f^z-probability one, poo is a well defined density matrix- 
valued random variable which does not depend on the value of pQ . Since the eigenvalue vector 
of /3o is hQ = (1,0, ...,0) G C'' we have that poo ~ t'feo- Now, the relation between the 
two families of measures is clear: the induced measure pd' is the distribution of the density 
matrix after one interaction, whereas z^tg is the distribution at the limit, after a large number 
of interactions. The reader may notice that this analogy is valid only in the case where 
h = (1,0, ...,0) (pure state on the environment). Generalizations of the (usual) induced 
measures to other environment states are possible, but out of the scope of the present work. 
To further compare the asymptotic and the one interaction induced measures, we plotted the 
spectra of samples of density matrices from both families in Figures [1] and [2j In particular, 
one should compare Figure 1(a) with Figure 2(a) \{d = d' = 2), Figure |l(d)| with Figure |2(b)| 



{d = d' = 3) and Figure 1(f) with Figure 2(c) {d = 3,d' = 5). 



^ure 


2(a) 


2(c: 


{d = 



One particularly simple case is obtained by taking b = {1/d' , . . . ,1/d'). The measure Vf, is 
then trivial, being equal to the Dirac mass supported on the "chaotic state" 1/d. In the next 
lemma we prove some basic properties of the newly introduced measures i^^. A more thorough 
investigation of these measures is postponed to a later work. 

Proposition 4.6. The probability measures Vf, have the following properties: 

(a) For every probability vector b, the measure i^b "is well defined, in the sense that the 
distribution of poo = linin^oo['&^''^]"(/Oo) does not depend on the eigenvectors of P, but 
only on the eigenvalue vector b. 

(b) For all unitary matrix V € hl{d), p and VpV* have the same distribution (we say that 
the measure Vf, is unitarily invariant^. 

(c) There exists a probability measure on the probability simplex Ad-i such that if D 
is a diagonal matrix sampled from Uf, and V is an independent Haar unitary on lA{d), 
then VDV* has distribution v^. In other words, the distribution of a random density 
matrix p ^ i^h is determined by the distribution of its eigenvalue vector Ad-i B X ^ ni^. 

Proof. To prove the first assertion, we show that for all W E lA{d'), replacing (3 with W f3W* 
does not change the distribution of poo- To see this, note that by the invariance of the 
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(d) (e) (f) 

Figure 1. Asymptotic measure eigenvalue distribution. First row, from left 
to right: {d = 2,b = [1,0]), (d = 2,b = [3/4,1/4]), {d = 2,b = [1,0,0,0]). 
Second row: (d = 3,6 = [1,0,0]), (d = 3,b = [3/4,1/8,1/8]) and (d = 3,b = 
[1,0,0,0,0]). 




eigenvalues 



(a) (b) (c) 

Figure 2. Induced measure eigenvalue distribution for (d = 2,d' = 2), (d = 
3,d' = 3) and {d = 3,d' = 5). 

Haar probability measure i)dd', the random matrices U and U = U{ld'S'W) have the same 
distribution. It follows that the same holds for the random channels and <I)^'^ = ^U,Wi3W* 
and thus for their invariant states. The second affirmation is proved in the same manner (this 
time using a fixed unitary V acting on Ti) and the third one is a trivial consequence of the 
second. □ 

5. Repeated interactions with random auxiliary states 

In the previous section we considered repeated identical quantum interactions of a system 
S with a chain of identical environment systems £. We now introduce classical randomness 
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in our model by considering random states on the environment £. In this model, the unitary 
describing the interaction is a fixed deterministic matrix U G U{dd'). 

The n-th interaction between the small system S and the environment E is given by the 
following relation: 

Pn = $^"(Pn-l) = Ttjc[U{pn-l ® Pn)U*], 

where (/?n)n is a sequence of independent identically distributed random density matrices. 
Notice that, since U is constant, we use the shorthand notation = 

We are interested, as usual, in the limit n — > oo. In this case however, the (random) 
channels do not have in general a common invariant state, so one has to look at ergodic 
limits. We use here the machinery developed by L. Bruneau, A. Joye and M. Merkli in [6] 
(see [71 H] for additional results in this direction). For the sake of completeness, let us state 
their main result. 

Theorem 5.1 ([6], Theorem 1.3.). Let {Mn)n be a sequence of i.i.d. random contractions of 
M.d{C) with the following properties: 

(a) There exists a constant vector ip € such that M[uj)ip = tp for (almost all) to; 

(b) F{the multiplicity of the eigenvalue 1 of M{lo) is exactly one) > 0. 

Then the (deterministic) matrix K[M] has eigenvalue 1 with multiplicity one and there exists 
a constant vector G such that 

1 ^ 

^ V Afi(c^)M2 (.;)•• -M, (a;) = \i^){9\ = Pi,e[m], 

N — >oc iV — 
n=l 

where -Pi,E[Af] the rank-one spectral projector o/E[M] corresponding to the eigenvalue 1. 

Note that this result does't apply to our situation, mainly for two reasons: the order of the 
composition of the channels ^ is reversed and the linear applications <I>^" do not necessarily 
share a constant invariant state Tp. This inconvenient can be overcome by considering dual 
channels (see Section [3]), or, in physicists' language, by switching from the Schrodinger to 
the Heisenberg picture of Quantum Mechanics. Duals of quantum channels are unital, hence 
they have in common the invariant element I. Another important benefit of considering duals 
is that the order of composition of maps is reversed. Indeed, if one starts from a state 
applies successively n channels ...,<!>„ and finally measures an observable A G A^^^(C), 
it is easy to see that the expected outcome is 

TV[($„ o . . . o $i)(p) . ^] = TY[($„_i o . . . o ^.i)(p) . m^{A)] = . . . = Tr[p • (^1 o . . . o *„)(^)]. 

We are now in position to state and prove the analogue of Theorem 15.11 for infinite products 
of quantum channels, simply by replacing quantum channels with their duals. 

Theorem 5.2. Let {^n)n be a sequence of i.i.d. random quantum channels acting on A4d{C) 
such that P(<I> has an unique invariant state) > 0. Then E[$] is a quantum channel with an 
unique invariant state 6 G A^^'^(C) and, W -almost surely, 

1 ^ 

lim - ^[<I>„ o . . . o ci>i](po) = 0, Vpo G MY{C). 

n=l 

Proof. Let us start by introducing some notation. Let, for some initial state po G AI^'^(C), 

1 ^ 

/^7V = ^ J^[^nO---0$l](po), 

n=l 
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and consider the dual operators which are, as described earlier, the adjoints of with 
respect to the Hilbert-Schmidt scalar product on M.d{C). Then, for a self-adjoint observable 
A e A^^'"(C), one has 



(7) TiIpnA] = Tr 



1 ^ 



n=l 



It is easy to see that the random operators ^„ satisfy the hypotheses of Theorem 15.11 on 
the Hilbert space A^rf(C) endowed with the Hilbert-Schmidt scalar product. Indeed, the 
spectrum of ^ is the complex conjugate of the spectrum of <I>, hence ^' is a contraction (with 
respect to the Hilbert-Schmidt norm). Moreover, with non-zero probability, 1^ is the unique 
invariant state of ^. From the Theorem 15. H one obtains the existence of a non-random 
element 9 G A4d{C) such that, P-almost surely, 

1 ^ 

hm ^vl/,o...ovl.„=|l,)(0|. 

n=l 

Plugging this into Eq. ([7]), one gets 

lim Tr[^^^] = Tt[po\U){9\A] = (0, A)hs Tr[po IJ = TT[e*A]. 

N~too 

Since the set of density matrices A^^'^(C) is (weakly) closed, 9 = 9* G ^A^'^{C) and 
lim^^oo f^N = 9. The fact that 9 is the unique invariant state of E[$] follows again from 
Theorem 15.11 □ 

Remark 5.3. When comparing the preceding theorem with the Proposition 14. 11 one notes that 
the hypotheses are relaxed here, asking only that the eigenvalue 1 is simple, without further 
constraints on the peripheral spectrum. This is due to the fact that we are considering Cesaro 
means and fluctuations (such as the ones seen in Example 14. 2p cancel out at the limit. 

We now move on to apply the preceding general result to the setting described in the 
beginning of this section. Recall that the successive interactions were described by i.i.d. 
random quantum channels = where 

^^{p)=TT!c[U{p®(3)U*]. 

Since the previous equation is linear in /3, E[$^] = ^'^I^l and one gets the following corollary. 

Corollary 5.4. Let {Pn}n be a sequence of i.i.d. random density matrices and consider the 
repeated quantum interaction scheme with constant interaction unitary U . Assume that, with 
non-zero probability, the induced quantum channel <I>^ has an unique invariant state. Then, 
F-almost surely, for all initial states po G A^^'^(C), one has 

J™ ^^[^^"o...o<I>/5i](po)=0, 

Af^oo TV 

n=l 

where 9 G A^^'^(C) is the unique invariant state of the deterministic channel ^'^I^l. In 
particular, i/E[/3] = 1^/ /d', then 9 is the "chaotic" state Id/d. 
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6. Repeated interactions with i.i.d. unitaries 

We now consider a rather different framework from the one studied in Sections H] and [H 
We shall assume that the interaction unitaries Un acting on the coupled system Ti ® IC are 
random independent and identically distributed (i.i.d.) according to the unique invariant 
(Haar) probability measure \]dd' on the group lA{dd'). This is a rather non-conventional model 
from a physical point of view, but it permits to relax hypothesis on the successive states of 
the environment and to obtain an ergodic-type convergence result. 

As before, we start with a fixed state po € 7W^'^(C). The n-th interaction is given by 
Pn = where (/3n)n is a (possibly random) sequence of density matrices on 

K and {Un)n is a sequence of i.i.d. Haar unitaries of U{dd') independent of the sequence 
{Pn)n- Note that we make no assumption on the joint distribution of the sequence {Pn)n] in 
particular, the environment states can be correlated or they can have non-identical probability 
distributions. The state of the system after n interactions is given by the forward iteration of 
the applications $'^'1'/^": 

(8) Pn = O $t^n-l:/3n-l q . . . O ^'^^''^^pQ. 

Since we made no assumption on the successive states of the environment /?„ G A1^;^(C), 
the sequence {pn)n is not a Markov chain in general. Indeed, the density matrices {Pn)n were 
not supposed independent, hence Pn+i (and thus Pn+i) may depend not only on the present 
randomness, but also on past randomness, such as (3n-i, Pn-2, etc. Although the sequence 
{pn)n lacks markovianity, it has the following important invariance property. 

Lemma 6.1. Let {Vn)n be a sequence of i.i.d. Haar unitaries independent of the family 
{Un,Pn}n cind Consider the sequence of successive states {pn)n defined in Eq. (0). Then the 
sequences {pn)n cind (VnPnV^)n have the same distribution. 

Proof. Consider a i.i.d. sequence (Ki)n of fi^-distributed unitaries independent from the Un's 
and the /?„'s appearing in Eq. ([8]). To simplify notation, we put p„ = VnPnVn- We also 
introduce the following sequence of (random) dd' x dd' unitary matrices: 

Ui = {Vi^l)Ui, 

Un = (K ® l)Un{V:-i ® I), Vn ^ 2. 
A simple calculation shows that 

p^ = ^U„,f3„ Q ^U„-i,p„-i o . . . o $t^i-/3ipQ. 

It follows that, in order to conclude, it suffices to show that the family (Un)n is i.i.d. and 
[)rfrf'-distributed (it is obviously independent of the /3's). We start by proving that, at fixed n, 
Un is f)drf'-distributed. Since the families {Un)n and (Vn)n are independent, one can consider 
realizations of these random variables on different probability space Un '■ ~^ hl{dd') and 
: f^n ^ U{d). For a positive measurable function / : U{dd') — > M+, one has (we put Vq = I) 

nf{Un)] = E[/((K ® l)Un{V:^l ® I))] = 

= / /((K(^^)®i)^n(u;i)(K:-i(^Li)®i))dn'^n)^n^n)rfn^n-i) 

= / (/ /((K(o;^)0l)^n(a;i)(K-l('^n-l)®I))^(^n)) dF{u;l)dF{ul_^) 

j nf{Un)W{u^i)dn^i.,) = nfiUn)], 
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where we used in (*) the fact that the Haar probabihty on lA{dd') is invariant by left and right 
multiphcation with constant unitaries. We now claim that the r.v. Un are independent. For 
some positive measurable functions /i, ...,/„ : U{dd') M+, one has 



E 



n ikiUk 



.k=l 



Ufkm^mkiVkU^^)) = 

k=l 

n n 

= I nA((^fe('^^)®i)f^fe('^^)(^fc-i(^Li)®i))n^('^^)'^^(^fe) 

k=l k=l 

= fflif fk{mu;l)®l)UkU){V,Uiul_,)®l))dF{u^^^^ 

k=i ^ k=\ 

„ n n n 



k=l 



k=l 



k=l 



Again, we used in the equality (**) the invariancc of the dd'-dimensional Haar measure and 
in (* * *) the fact that and 17^ have the same distribution. 

□ 

We conclude from the above result that although the successive states of the small system 
{Pn)n are random density matrices that can be correlated in a very general way, their joint 
probability distribution is invariant by independent unitary basis changes. In other words, 
the correlations manifest only at the level of the spectrum, the matrices being independently 
rotated by random Haar unitaries. The ergodic convergence result in such a case is established 
in the following proposition. 

Proposition 6.2. Let (Tn)n be a sequence of random density matrices (we make no as- 
sumption whatsoever on their distribution) and {Vn)n o- sequence of i.i.d. Haar unitaries 
independent of (r„)„. Then, almost surely, 



VmV,* + ... + VnTnV* 



n 



d' 



Proof. Since both sides of the previous equation are self-adjoint matrices, it suffices to show 
that for any self-adjoint operator A £ 7Wrf(C) we have lim„^oo Tr[cr„yl] = Tr[A]/(i. Using 
the invariance of the Haar measure, one can assume that the observable A is diagonal A = 
Yli=i ™ some fixed orthonormal basis {ei}f^^ of C*^. In the same basis, we write 

Tk = (t-j and Vfe = {v^j)fj^i. To simplify notation, we put 

n 

where n = Tr[VkPkV*A] = Et^,j=i Cl^^^S^S" Using the fact that 



E 



(k) (fe) 



one can easily check that that the random variables have mean Tr[^]/d, finite variance 
(a rough bound for E[r|] is Tr[^]^) and that they are not correlated (cov(rj., T^./) = 0, if 
k ^ k'). It is a classical result in probability theory that in this case the (strong) Law of 
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Large Numbers holds and thus, almost surely, 

lim TrianA] = 

n—KX a 

□ 

Putting the previous proposition and Lemma 16.11 together, one obtains the main result 
of this section, an ergodic-mean convergence result for the sequence of states of the "small" 
system. 

Proposition 6.3. Let {pn)n the successive states of a repeated quantum interaction scheme 
with i.i.d. random unitary interactions. Then, almost surely, 

pi + ...+ Pn Id 



lim 



n d 
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